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J> ■ Abstract 

■ A global existence theorem, with respect to a geometrically defined 

| time, is shown for Gowdy symmetric globally hyperbolic solutions of 

£N| ■ the Einstein- Vlasov system for arbitrary (in size) initial data. The 

spacetimes being studied contain both matter and gravitational waves. 

OO ' 
OS 

O '. 1 Introduction 

cr 

*h . An important problem in classical general relativity is the question of global 

existence (in an appropriate sense) for globally hyperbolic solutions of the 
vacuum-Einstein and matter-Einstein equations. The main motivation being 
^ . its relationship to the cosmic censorship conjectures. Strong cosmic censor- 

ship has eg. by Eardley and Moncrief [EM] been formulated as a question 
on global existence and asymptotic behaviour of solutions to the Einstein 
equations, suggesting a definite method of analytical attack. 

To begin studying the long-time behaviour of solutions to a complicated 
partial differential equation system one might focus on families of solutions 
with some prescribed symmetry. With the exception of the monumental 
work on global nonlinear stability of the Minkowski space by Christodoulou 
and Klainerman [CK] , the practice in general relativity has for long been to 
study "global existence" problems under symmetric assumptions. 

One family of (cosmological) solutions which have been studied exten- 
sively is the Gowdy spacetimes [G] . These spacetimes are vacuum but admit 
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gravitational waves (in contrast to eg. spherically symmetric spacetimes). 
Global existence has been shown for the Gowdy spacetimes [M], strong cos- 
mic censorship is settled in the case of polarized Gowdy spacetimes [CIM], 
and much is known about the subset of the Gowdy spacetimes which admit 
an extension across a Cauchy horizon [CI]. 

In this paper we show global existence, with respect to a geometrically 
defined time, for matter spacetimes (Einstein- Vlasov) with Gowdy symme- 
try and thereby we extend Moncrief's result [M] in the vacuum case. This 
is the first result which provides a global foliation of a spacetime containing 
both matter and gravitational waves. Moreover, for matter spacetimes there 
are only a few global results available all together. Let us briefly mention 
some of these results. First, by matter spacetimes we have in mind space- 
times where the matter consists of massive particles. One can also consider 
spacetimes which only contains radiation and important results have been 
obtained in this direction, e.g. Christodoulou has obtained strong results 
in the spherically symmetric case with a scalar field as matter model (see 
e.g. [Cul], [Cu2] and the references therein). For spacetimes containing 
massive particles the main global results can be summarized as follows. Un- 
der a smallness condition on the initial data, Rein and Rendall have [RR] 
shown that solutions of the spherically symmetric Einstein- Vlasov system 
are geodesically complete. Some information on the large data problem 
was then obtained in [RRS]. Christodoulou has in a series of papers (see 
[Cu3] and the references therein) studied the Einstein-Euler equation in the 
spherically symmetric case for a special equation of state, adapted to un- 
derstand the dynamics of a supernova explosion. He can globally control 
the solutions to the Cauchy problem and he finds solutions whose behaviour 
resembles qualitatively that of a supernova explosion. Finally, the most rel- 
evant results in the context of this paper are those on cosmological solutions 
by Rendall [Rll-2] and Rein [Rn]. These are discussed in some detail in 
relation to our result below. 

Our method of proof is inspired by a recent global foliation result for 
vacuum spacetimes admitting a T 2 isometry group, acting on T 3 spacelike 
surfaces [BCIM]. These spacetimes are more general than the Gowdy space- 
times: both families admit two commuting Killing vectors but in the Gowdy 
case there is the additional condition that the twists are zero. The twists 
are defined by 

ci = e^sX^VPX 5 , c 2 = e^X^Y^Y 5 , (1) 
where X, Y are Killing vectors associated with the isometry group. It follows 
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from the Einstein equations that in vacuum these quantities are constant 
throughout spacetime [G]. 

One difficulty in studying long-time existence problems in general rela- 
tivity is the lack of having a fixed time measure. A solution which remains 
regular for an infinite range of one time scale may become singular within 
a finite range of another. In [BCIM] this problem is treated by choosing a 
coordinate system in which the time is fixed to the geometry of spacetime. 
In fact, the time is defined to be the area of the two dimensional spacelike 
orbits of the T 2 isometry group. These coordinates are called areal coordi- 
nates. The main theorem in [BCIM] shows that the entire maximal globally 
hyperbolic development of the initial hypersurface can be foliated by areal 
coordinates. These coordinates are however only used in a direct way in the 
future direction. To show that the past of the initial hypersurface is cov- 
ered by areal coordinates the authors use conformal coordinates (the time is 
not fixed to the geometry of spacetime) in which the equations take a more 
suitable form for an analytical treatment. By a long chain of geometrical 
arguments it is then shown that the development in conformal coordinates 
admits a foliation by areal coordinates, and that it covers the past maximal 
globally hyperbolic development of the initial hypersurface. 

We prove that T 3 x R— matter spacetimes with Gowdy symmetry admit 
global foliations by areal coordinates. The matter content is described by 
the Vlasov equation. This is a kinetic equation and gives a statistical de- 
scription of a collection of collisionless "particles" . In the cosmological case 
the particles are galaxies or clusters of galaxies whereas in stellar dynamics 
they are stars. The Vlasov equation has shown to be suitable in general 
relativity for the study of the long-time behaviour of matter in gravitational 
fields. In particular it rules out the formation of shell-crossing singularities. 
For a discussion on the choice of matter model see [R14] and [R15]. 

To prove the existence of a global folitaion we also work directly in areal 
coordinates in the expanding (future) direction, and in the contracting (past) 
direction, we first show a global existence theorem in conformal coordinates 
and then we invoke the geometrical arguments in [BCIM] to complete the 
proof. We point out that our result depends strongly on the exact strucure 
of the Vlasov equation and do not hold for general matter models which 
are only restricted by certain inequalities on the components of the energy- 
momentum tensor. 

A related and interesting result has recently been shown by Rendall [Rll] 
(see also [R12]). He considers T 2 symmetric spacetimes for the Einstein- 
Vlasov and the Einstein-wave map equations and he shows that if such a 
spacetime admits at least one compact constant mean curvature (CMC) hy- 
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persurface then the past of that surface can be covered by a foliation of 
compact CMC hyper surf aces. The CMC- and the areal coordinate folia- 
tion are both geometrically based time foliations which provide frameworks 
for studying strong cosmic censorship and other global issues. The main 
motivation for developing techniques to obtain CMC foliations is that the 
definition of a CMC hypersurface does not depend on any symmetry as- 
sumptions and it is hence possible that CMC foliations will exist for rather 
general spacetimes. The areal coordinate foliation used here is less general 
since it is adapted to the symmetry, but leads in the Gowdy case (note that 
the results in [Rll] apply to the more general T 2 symmetric spacetimes which 
we hope to explore in the future) to stronger results. Namely, the arguments 
in [Rll] do not show that the entire future of the initial hypersurface can be 
covered, and the existence of the CMC foliation is only guaranteed under 
the hypothesis that spacetime admits at least one such hypersurface. 

We also mention a result in this direction due to Rein [Rn]. He has 
studied cosmological Einstein- Vlasov spacetimes with stronger symmetry 
restrictions than in the Gowdy case (the spacetimes admit three Killing 
vectors). In these spacetimes gravitational waves cannot exist. For plane 
symmetry (the relevant case for us) he has shown existence back to the 
initial singularity for small initial data, and under the assumption that one 
of the field components is bounded, he obtains global existence for large data 
in the future direction. An interesting result in his work is that the initial 
singularity is shown to be a curvature singularity as well as a "crushing" 
singularity (see [ES]). 

The outline of the paper follows in the large that of [BCIM]. In section 
2 we describe Gowdy symmetry and give the equations for the Einstein- 
Vlasov system in areal and conformal coordinates. The main theorem is 
formulated in section 3 where we also describe the geometrical arguments in 
[BCIM] needed to complete the proof in the contracting direction. Section 
4 is devoted to the analysis in the contracting direction. Estimates for the 
field components and the matter terms are derived in conformal coordinates, 
by using e.g. light-cone arguments and methods originally developed for the 
Vlasov-Maxwell equation. The analysis in the expanding direction is carried 
out in areal coordinates in section 5 where a number of estimates are derived. 
Light-cone arguments and an "energy" monotonicity lemma are important 
tools for obtaining bounds on the field components and their derivatives. 
The control of the matter terms and their derivatives rely on three lemmas. 
The first one is the "energy" monotonicity lemma just mentioned. Then, in 
the second lemma a careful analysis of the characteristic system associated 
with the Vlasov equation is carried out, which leads to a bound on the sup- 
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port of the momenta. The third lemma provides bounds on the derivatives 
of the matter terms and relies indirectly on the geodesic deviation equation. 
This equation relates the curvature tensor and the acceleration of nearby 
geodesies and has proved useful in previous studies of the Einstein- Vlasov 
system (see [RR], [Rn] and [R13]). 

2 The Einstein- Vlasov system with Gowdy sym- 
metry 

Let us begin with a brief review of Gowdy symmetry. Consider a space- 
time that can be foliated by a family of compact, connected, and orientable 
hypersurfaces. If the maximal isometry group of the spacetime is two di- 
mensional, and if it acts invariantly and effectively on the foliation, then the 
isometry group must be U(l) x U(l). Moreover, the foliation surfaces must 
be homeomorphic to T 3 , S 1 x S 2 , S s or L(p, q) (the Lens space), and the ac- 
tion is unique up to equivalence. The Killing vector fields X, Y associated 
with the isometry group have to commute in such a spacetime. We say that 
spacetimes satisfying the symmetry conditions above and in which both the 
twists ci , C2 (see |l|) vanish have Gowdy symmetry. We remark that the term 
"Gowdy spacetime" is reserved for the vacuum case. For more background 
on Gowdy symmetry we refer to [G], [CI]. 

As mentioned above there are several choices of spacetime manifolds 
compatible with Gowdy symmetry. In this paper we restrict our attention 
to the T 3 — case. It is an interesting fact that in vacuum this is the only 
possibility if the condition of vanishing twists is relaxed. 

The dynamics of the matter is governed by the Vlasov equation. This 
is a kinetic equation and models a collisionless system of particles, i.e. the 
particles follow the geodesies of spacetime. For a nice introduction to the 
Einstein- Vlasov system see [R13] . We also mention the survey of Ehlers [E] 
for more information on kinetic theory in general relativity, and the book by 
Binney and Tremaine [BT] for some applications of kinetic theory in stellar 
dynamics. 

We will use two choices of coordinates, areal coordinates and conformal 
coordinates. It has been shown in [CI] that, at least locally, any globally 
hyperbolic (non-flat) Gowdy spacetime on T 3 x 1 admits each of these co- 
ordinates. Both sets of coordinates are chosen so that 
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are Killing vector fields (a, b, c and d are constants with ad — be / 0), and in 
both cases 9 <E S 1 denotes the remaining spatial coordinate. Below the form 
of the metric and the Einstein- Vlasov system is given in areal and confor- 
mal coordinates. The functions R, a, U, A, r/ all depend on t and 9 and the 
function / depends on t, 9 and »el 3 . 



Areal Coordinates 

Metric 

g = _ e 2(v-U) adt 2 + e 2(v-U) d9 2 + e 2U {dx + My) 2 + (2) 

The Einstein-matter constraint equations 

f = U? + aUl + ^{Al + aAl)+^- u ) ap (3) 
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= 2ta 2 e 2 ^- u) (P 1 - p) (5) 



The Einstein-matter evolution equations 
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mi-owe = — + ^-^ + — - U t+ aU e + ^2-( A t- aA e) 

- ae Kv-u)p 3 _ ^ ae 2{, + u)p 2 _ ?A ae 2v S23 (6) 

TT 11 Ut i Udaf) , Utat i e4f/ //l2 /l2N 

Efe-aEfo = - T + _ + _ + _(^-a^) 

+ I e 2(^) a(/0 _p 1+ p 2 _p 3) (7) 
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The Vlasov equation 



df yWd/ 
dt v° 38 



+^({v 3 ) 2 - {v 2 ) 2 ) 



U t V 2 + y/aU g 



v 1 v 2 



1 



(- - U t )v 3 - ^U ^ r + 



v 

v l y 3 e 2 U v 2 v 1 







i 2u 
T 

df 



(r] e -Ug + —)Vav +{rit-Ut)v 

df_ 

dv 1 



dv 2 

= (9) 

dv 3 1 ' 



The matter quantities 

P (t,o) = [ v°f(t,e,v)dv (io) 

Pk(t,0) = [ ^Q-f(t,8,v)dv, k = 1,2,3 (11) 
J(t,0) = f v^f&O^dv (12) 

/* 2 3 

S23M) = / ^-/(t,f,«)d« (13) 

Here the variables v are related to the canonical momenta p through 

v o = ^ e v- V, t; 1 = efo-V, « 2 = eV + Ae u p 3 , v 3 = te~ V, (14) 
and 

p»:=^, x» = {t,0,x,y), 

where r is proper time. It is assumed that all 'particles' have the same 
mass (normalized to one) and follow the geodesies of spacetime (collisionless 
particle system). Hence 

Wp" = -i, 

so that 

v° = v / l + (^ 1 ) 2 + (« 2 ) 2 + (v 3 ) 2 . (15) 

In conformal coordinates the function a is removed, having the conse- 
quence that the orbital area function R now depends on both t and 8 (in 
areal coordinates R = t). In these coordinates the metric and the Einstein- 
Vlasov system take the following form. 
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Conformal coordinates 



Metric 



g = e 2(v-U) { _ dt 2 + + e 2U {dx + My) 2 + e -2U R 2 dy 2 



(16) 



The Einstein-matter constraint equations 



u + ui + 



^ 2 + ^ + T 



Roe rj t R t rjeRe 



R 



R 



,2{n-u) 



P (17) 



2U t U e + e — 2 A t A e + ^ - ^ - ^ = e 2 (' n ~ u ^J (18) 

(19) 



The Einstein-matter evolution equations 



^tt - ^ 
Ptt — Pee 
»7tt - Vee 



U e R e U t R t , e 



4C7 



P P 2P 2 

+ I e 2(^) (p _p i+ p 2 _p 3 ) 



2 

P*^ 



(20) 



^ + 4(^C/ e - ^) + 2Re 2 ^U)s 23 (21) 

-ft XI 



Pe^^-P^ 



Ul-U? + ^- 2 {Al-Al)-^-v) P , 
_ge 2 ^)p 2 -^e 2 ^ 23 



(22) 



(23) 



The Vlasov equation 



dt v° 86 



( Ve - U e )v° + ( Vt - U t )v l - U e ^f- 



2, ,31 



, (TT ReAv s Y A e 2U v~v 



(f - U t )v 3 - (U e 



1„,3 



Rg V l V 
~R>~ 



dv 1 

e 2U y 2 



U t v 2 + U e - 



v 1 v 2 ~ i 



dv 2 



df_ 

dv 3 



(24) 
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The matter quantities p,Pk,J and S23 are given by (|l0|)-(|l3|), where in this 
case 



v 



° = e "- V, v 1 = efo-V, v 2 = eV + Ae u p 3 , v 3 = Re~ u p 3 , (25) 



and dig ) holds here as well. 

Remark. It might be instructive to relate the metric in (16) with that 
used by Rein [Rn] mentioned in the introduction. By letting A = and 
U = (1/2) In R in (|l6|) we obtain a metric which admits three Killing vectors 
and which depends on two field components. The distribution function / 
depends in this case on p 1 and (p 2 ) 2 + (p 3 ) 2 - 



3 The main theorem 

Let (h, k, fo) be a Gowdy symmetric initial data set on T 3 , by this we mean 
that h is a Riemannian metric on T 3 , invariant under an effective T 2 action; 
k is a symmetric 2-tensor on T 3 , also invariant under the same T 2 group 
action; the twists c\ and C2 are both zero; the initial distribution function fo 
is defined on T 3 and is invariant under the same T 2 group action and pos- 
sesses the following additional symmetry, which reads, in coordinates that 
cast the metric in the forms (|2|) or (|ilf), fo{6,p l ,p 2 ,p 3 ) = fo(&,P l ,—p 2 ,P 3 ) = 
foiOiP 1 ,p 2 , — p 3 ) (this assumption is necassary for the Einstein- Vlasov sys- 
tem to be compatible with the form of the metric); and that (h, k, fo) satisfy 
the Einstein- Vlasov constraint equations. We also assume that (h, k) are C°° 
on T 3 and that fo is a non-negative, not identically zero, C°° function of 
compact support on the tangent bundle T(T 3 ) of T 3 . 

Remark. The smoothness assumption on the initial data is not a necas- 
sary condition. It is included so that we can refer directly to the classical 
local existence theorems. However, the estimates in this paper provide the 
information needed for proving a local existence theorem for C 2 x C 1 data 
(h, k) and C 1 data fo. Moreover, the assumption fo ^ is here included 
for a technical reason and we refer to [M] or section 5 in this paper for the 
vacuum case. Indeed, it is in this case possible to work directly in areal 
coordinates and the estimates derived in section 5 are sufficient. See also 
the remark following lemma 1 in that section. 

The results by Choquet-Bruhat [CB] and Choquet-Bruhat and Geroch 
[CBG], show that there exists a unique maximal globally hyperbolic de- 
velopment (S x R, g, f) of a given initial data set on a three-dimensional 
manifold S for the Einstein- Vlasov equation. Let us briefly comment upon 
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the initial conditions imposed. The relations between a given initial data set 
(h, k) on a three-dimensional manifold £ and the metric g on the spacetime 
manifold is that there exists an imbedding if; of £ into the spacetime such 
that the induced metric and second fundamental form of ip(T,) coincide with 
the result of transporting (h, k) with ip. For the relation of the distribution 
functions / and /o we have to note that / is defined on the mass shell (for 
m = 1 it is the set of all future pointing unit timelike vectors). The initial 
condition imposed is that the restriction of / to the part of the mass shell 
over should be equal to f^o , d{il))~ l ) o <f> where </> sends each point 
of the mass shell over ip(T,), to its orthogonal projection onto the tangent 
space to VK^)- 

Our main theorem now reads, 

Theorem 1 Let (h,k,fo) be a smooth Gowdy symmetric initial data set 
on T 3 . For some non-negative constant c, there exists a globally hyperbolic 
spacetime (A4,g,f) such that 

(i) M = (c,oo) x T 3 

(ii) g and f satisfy the Einstein- Vlasov equation 

(in) M. is covered by areal coordinates (t,8,x,y), with t £ (c,oo), so the 
metric globally takes the form (^). 

(iv) (M,g,f) is isometrically diffeomorphic to the maximal globally hyper- 
bolic development of the initial data (h,k,fo). 

As described in the introduction we prove global existence in conformal 
and areal coordinates for the past and future directions respectively. Then, 
in order to prove Theorem 1 in the past direction, we need to invoke sub- 
stantial geometrical arguments from [BCIM]. For the future direction only 
a simple geometrical argument is needed for completing the proof. It should 
be pointed out that even if the geometrical results in [BCIM] concern the 
vacuum case they are true also for matter spacetimes as long as the Einstein- 
matter equations form a well-posed hyperbolic system, which of course is the 
case here. 

In section 4 we show that the past maximal development of (h, k, /o) 
in terms of conformal coordinates, which we denote by D~ on j:(h, k, /o), has 
t — > — oo as long as R stays bounded away from zero. Starting from this 
result we briefly describe how the geometrical arguments in [BCIM] lead to 
a proof of Theorem 1 in the past direction. 

First, in [BCIM] R is shown to be positive everywhere in the globally 
hyperbolic region of a T 2 symmetric spacetime. Also, along any past in- 
extendible timelike path in D~ on j(h, k, /o), R is shown to approach a limit 
Ro ^ (to be identified with c in Theorem 1), which is independent of the 
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choice of path. Moreover, for any R G (Rq,Ri), where i?i is the minimum 
value of R on the initial hypersurface, the level set R = R in D~ on f(h, k, fo) 
is shown to be a Cauchy surface. From these facts it follows from argu- 
ments in [CI] that D~ on f(h, k, fo) admits areal coordinates to the past of 
the hypersurface R = R\. Proposition 4 and 5 in [BCIM] then show that 
D~ on j(h, k, /q) is also isometrically diffeomorphic to the maximal past de- 
velopment, D~ (h, k, /o) of (h,k,fo) on T 3 . 

In the future direction, global existence in areal coordinates is almost 
sufficient for proving Theorem 1. The only statement that remains to be 
proved in Theorem 1 is that the future maximal development is covered by 
Ah, k, fo). This follows from a very short geometrical argument given 
in the proof of Proposition 5 in [BCIM]. 

4 Analysis in the contracting direction 

The local existence theorem of Choquet-Bruhat [CB] together with the result 
of Chrusciel (lemma 4.2 in [CI]) imply that for any Gowdy symmetric initial 
data set (h,k,fo) on T 3 , we can find an interval (ii,^) and C°° functions 
R,U,r] on (ii,i 2 ) x T 3 , and a non-negative C°° function / on (ii,^) x P 
(P denotes the mass shell) such that: these functions satisfy the Einstein- 
Vlasov equations in conformal coordinate form and for some to G (£1^2)5 
the metric g induces initial data on the to-hypersurface which is smoothly 
spatially diffeomorphic to (h,k), and the relation between / and fo given 
above holds. 

Now, in order to show that D~ on j(h,k, fo) has t — ► —00, as long as R 
stays bounded away from zero, it is sufficient to prove that on any finite 
time interval (t, to], the functions R,U, A,rf, f and all their derivatives are 
uniformly bounded and that the supremum of the support of momenta at 
time t, 

Q{t) := sup{M : 3(s, 6) G [t, t ] x SSuch that/(s, 6, v) / 0}, (26) 

is uniformly bounded. Note that the last condition implies that the matter 
quantities and their derivatives are uniformly bounded (if \df /dx^\ < C). 

Step 1. (Monotonicity of R and bounds on its first derivatives.) 
This is a key step and relies on Theorem 4.1 in [CI] together with the ar- 
guments in [BCIM]. We have to check that the matter terms have the right 
signs so that these arguments still hold. The bounds on R and its first 
derivatives will play a crucial role when we control the matter terms below. 
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First we show that \7R is timelike. Let us introduce the null vector fields 
dz = -L(d t + d e ), d x = ^=(d t -d e) , (27) 

and let us set = d^F, F\ = d\F for a function F. After some algebra it 



follows that the constraint equations (17) and (|lq) can be written 

deRi = rjzRt - RUj 1 - Re^-U) ( p -J), (28) 
d e R x = i lx R\ ~ RUl - Re 2 ^- U \p + J). (29) 

Let hi and h 2 be defined by 

h x ■= RU 2 + Re 2 ^~ u \p - J), 

and 

h 2 :=RUi + Re 2 ^- u \p+J). 
From ( JTOj ) and ( |T2| ) we have p > \J\, and since R > it follows that both 



hi and h 2 are non-negative. Solving equation (|28|) gives for any 9q E [0, 2tt] 
(suppressing the i-dependence) 

R^e) = A v ^ a)da R^e ) - f ^ ^ a)da hi{e)de. (30) 

Since R is C°° on S 1 it can be identified with a periodic function on the 
real line. If now R^(9q) = for any #0 then R^(2tt + 60) = 0, but from (|30|) 
this is only possible if hi vanishes identically. However, in the non-vacuum 
case (recall fo 7^ 0) h^(t,-) is strictly positive on some open set of [0, 2-k]. 
Therefore R^ is nonzero and has a definite sign. The same arguments apply 
to R\, and it follows that 

g^d^Rd p R = e-^-^R^Rx 

is strictly positive or strictly negative. The former possibility is ruled out 
since dgR = at some point on S . Thus \7R is timelike. This means that 
dtR is nonzero everywhere. Our choice of time corresponds to contracting 
T 2 orbits so that d t R > 0. 

Next we show that dtR and \dgR\ are bounded into the past. The evo- 
lution equation ( |22|) can be written 

d x Rt = Re 2 ^ u \p-Pi), (31) 
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or equivalently, 

d (i Rx = Re 2(r >- U) (p-Pi)- (32) 



The right hand side is positive since p > P\, see (10) and (11), and from 



(31) it follows that if we start at any point (to, 9q) on the initial surface we 
obtain 

R^e + s,t -s) <Rt(t ,6 ), (33) 
and similarly from (|32]) 

Rx(9 -s,t -s)<R x (t ,e ). (34) 

From these relations we get for any t £ (t, to) and any 6 £ S 1 , 

Rt(t,6) < sup R^(t ,e), (35) 
ees 1 

R x (t,9) < sup R x (t ,9). (36) 
ess 1 

This yields 

Rt(t,6) < swp(Rt + R x )(t ,e), (37) 

6&S 1 

and since VR is timelike everywhere we have \Rt\ > \Re\ and we find that 
both R t and \Rg\ are bounded into the past, so R is uniformly C 1 bounded 
to the past of the initial surface. 

Step 2. (Bounds on U, A and rj and their first derivatives.) 
The bounds on Ut,A t ,Ug and Aq to the past of the initial surface are ob- 
tained by a light-cone estimate, which in this case, with one spatial dimen- 
sion, is an application of the Gronwall method on two independent null 
paths. Then, by combining these results, one obtains the desired estimate. 

The functions involved in the light-cone argument are quadratic func- 
tions in the first order derivatives of U and A, defined by 

1 p 4t/ 

G = -R(U? + U e 2 ) + —(A 2 t +Al), (38) 

H = RU t U e + —A t A e . (39) 

We will see below that if we let the vector fields along the null paths act on 
G + H and G — H we obtain equations appropriate for applying a Gronwall 
argument. First, however, we motivate the definition of G and H by showing 
that they are components of an "enegy-momentum tensor" of a wave map. 
So let us follow [BCIM] and give a brief discussion of such a map in this con- 
text. Consider a base Lorentzian manifold (JV, v) with the two dimensional 
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manifold J\f corrseponding to the past conformal coordinate development of 
(h, k, fo) with the metric 



v := -dt 2 + de 2 , 



and consider a family of target manifolds (M. 2 ,d( t ^) with 



d {tfi) : =R(t,0)dU 2 + ^ F) dA 2 . 

Define the map 

$ : M — > R 2 , 

(t,6) ^ $(t,6) = (U(t,9),A(t,6)), (40) 

with du g\ providing an inner product on their tangents, e.g. 

< $ t , $ >= d ab ^e = RUtUe + —A t A g . 

There is a covariant derivative D compatible with u and semi-compatible 
with d( t fi) . Using Greek indices for the base and Latin indices for the target 
the action D can be expressed as 

d^i = d^i + - r^K, (4i) 

with 



^UU — °> ^UA — °> ^AA 



c 



4U 



2R 2 ' 



^AA — 0) ^AU — 2 > ^UU — 0) 

and the base Christoffel coefficients = (in the expanding direction 
these will be nonzero). As pointed out above, D is compatible with the flat 
metric v but not with d^g), we have 

p 4(7 

D x d ab = i? A (« - W)- 

By defining the wave operator □ := v^^D\D^ the evolution equations for 
U and A take the form 



UA = -^i + ^i + 2 Re 2 ^S 2 



R R ' ~~ 
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Now we define an "energy-momentum tensor" for the maps <£, 

AU i 4C7 

= RU X U, + ^ A A A„ + -^A M [i?(f/ f 2 - U$) + —(A 2 - A 2 e )}. 

We find that G = T tt and H = T te . 

Let us now derive bounds on U and A and their first order derivatives. 
By using the evolution equation (|20|) and (|2l|) we find 



1 / p w 
-^-rJu?-U 2 + ^(-A 2 ' 



d x (G + H) = —=R* Uf-U e z + —,(-Ai + A 



+ ^ U ^(v-u) (p _p 1 + p 2 _p 3] 



and 



-1 



3 4[/ 



+ K Uxe 2(V-U)( p _p 1+ p 2 _p 3) 

Now, integrating these equations along null paths starting at (t\,9) and 
ending at the initial — surface, and adding the results we obtain 

G(h,6) = \[G + H](t ,6 - (t - ti)) + \[G + H](t , 9+ (to- h)) 

'to 



-\ 1° K 1 (s,9-(s-t 1 )) + K 2 (s, 6 + (s- h)) ds 
1 Jti 

l - J ° [U^T] (s,9-( S -h)) + [U X T] ( s ,e + (s-h)) ds 
2 J [^^T](s,e-(s-t 1 )) + [—A x T]( S ,0 + (s-t 1 ))ds, (42) 



2 

where we have introduced the notations 



-1 / p 4U 
Ki = ^R X [(U 2 -U 2 + -jpirA + A 2 e ) ) , ( !.">) 
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-1 / P 4U \ 

K 2 = —R i i^V 2 - V 2 + —{-A* + A 2 )) , (44) 

T = ^- u Hp-P 1 + P 2 -P 3 ) (45) 
f = Re 2 ^) S 23 . (46) 

Let us first consider the matter terms. Note that for any t £ (t,to), the 
evolution equations ( |3l[ ) and ( |32| ) give 

i?g(to, 9 + (t - t)) - R^t, 9) = V2 J'" [Re 2 ^-^ (p - Pi)](s, 9 + (s - t))ds, 

(47) 

and 

R x (to,0-(t -t))-Rt(t,e) = V2 [ t0 {Re 2 (v-u\ p _ Pl)] ( s ^_ {s _ t))ds . (48) 



Hence, since i? is uniformly C 1 bounded to the past of the initial surface it 
follows that the right hand sides are uniformly bounded. From QT0|)-(11) we 

have 

P> P 1 +P 2 + P 3 , 

and thus 

< (p - Pi + P 2 - P 3 ) < 2{p - P x ), 

and from (|0|) and the elementary inequatlity 2ab < a 2 + b 2 , a, b G M, we 
have 

2|5 23 | <P 2 + P3<P-Pi. 
In view of ( |47| ) and (|4^) we therefore have that both 

J° T(s,0±{s-t))ds, (49) 

and ^ 

J° \f(s,e±(s-t))\ds, (50) 

are uniformly bounded on (t,io] x sa Y by C(i). Now, by using the 
inequality 2ab < a 2 + b 2 again, we get 



2G\ 
~R) 



1/2 
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and 

e 2U , , , /2C^ 1/2 



2R l 51 ~ V R 

The same estimates also hold for U\ and A\. Since and R\ are uniformly 
bounded it clearly follows that 

■ . CG i , GG 

for some constant G. So the identity (42) now implies, as long as R stays 
bounded away from zero, that 

supG(h,e) < sup2G(t ,#) +sup2G(t ,6) (51) 
e e e 

)ds. (52) 



+G sup (G(-,-)) 1/2 + /°G sup G(s, 

[hMxS 1 Jti 9 



On a sub time interval [cxi, Co] C [ii, io] on which sup s i G(s, •) is increasing 
as s — > ct"i, which is the case of interest for us since we wish to obtain a 
uniform upper bound on G, we have 

sup (G(-,-)) 1/2 < S up(G(a 1 ,-)) 1/2 . 

[o-i.o-olxS 1 S 1 

Therefore, by dividing both sides by supgi (G(a\, •)) 1 ^ 2 we obtain a Gronwall 
inequality on the time interval [<ti,<to] for sup^i G 1//2 leading to a upper 
bound on sup s i G(a±,-) as long as R stays bounded away from zero. By 
repeating this argument we conclude that sup 5 iG is uniformly bounded on 
(t,to\, leading to bounds on U and its first order derivatives, and thus also 
on A and its first order derivatives, as long as R stays bounded away from 
zero. The bounds on |?7|, |rft| and 1 77^ | are obtained in a similar way since the 
evolution equation (p3h can be written 



d xV , = ^-L/ 2 + |^(^ 2 -^)-e 2 ^^(P3 + ^P 2 + ^e 2 ^ 23 ), (53) 
or equivalently, 

%A = ^-^ 2 + |^(^-A 2 )-e 2 (''-^(P3 + ^P 2 + Me 2 ^ 2 3). (54) 

We found above that the integrals along null paths for the matter quantity 
— Pi) were bounded to the past of the initial surface. Therefore, 
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since < Pk < p — Pi, k = 2, 3 and IS23I < p — Pi we have, as long as R 
stays bounded away from zero, that the integrals along the null paths for 
the matter terms in the rigt hand sides above are bounded as well, since U 
and A are bounded. Now, since the first order derivatives of U and A are 
uniformly bounded we immediately obtain that \rj^\ and |^| are bounded 
by integrating the equations for r\ along null paths. Since r\ t = 775 (?7f + ?7a) 
and rjt = ^= (rj£ — rj\) we find that 77 is uniformly C 1 bounded to the past of 
the initial surface as long as R stays bounded away from zero. 

Step 3. (Bound on the support of the momentum.) 
Note that a solution / to the Vlasov equation is given by 

f(t, 9, v) = / (6(0, t, 9, v), V(0, t, 6, v)), (55) 

where and V are solutions to the characteristic system 

d& _ V^_ 

dV 1 (V 2 \ 2 
^_ = - {m - Ud )v^-{ m -U t )V l + U e K -^l- 

Re (V 3 ) 2 A e m V 2 V 3 
~ {Ue ~ -R ] — + li e -yo"' 

dV 2 „„, „v l v 2 



U t V< - Ug- 



ds 1 ° V° ' 

dV 3 Rt uW 3, (n Re^V 3 

e 2U , . , V 1 



-^(At + Ae^V 2 , 

and Q(s,t,x,v), V(s,t, x,v) is the solution that goes through the point 
(6, v) at time t. Let us recall the definition of 

Q(t) := sup{\v\ : 3(s,0) G [t,t ] x ^such that/(s, 9, v) / 0}. 



If Q(t) can be controlled we obtain immediately from (|10|)-(ll2|) bounds on 
p, J, 5*23 and P k , k = 1,2,3, since \\f\\oo < 11/olU from (g). Now, all of 
the field components and their first derivatives are known to be bounded on 
it, to]; a s long as R stays bounded away from zero. Also, the distribution 
function has compact support on the initial surface and therefore |y fc (to)| < 
C. So by observing that \V k \ < V°, k = 1,2, 3, a simple Gronwall argument 
applied to the characteristic system gives uniform bounds on |y fc (t)|, t £ 
(i,io]> an d it follows that Q(t) is uniformly bounded on (t, to]. 
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Remark. By a Killing vector argument, bounds on \V 2 \ and | X^ 3 1 can 
be derived if merely \U\ and are bounded and R > e > 0. Such an 
argument will be used in the expanding direction. 

Step 4- (Bounds on the second order derivatives of the field components 
and on the first order derivatives of /.) From the Einstein- matter constraint 
equations in conformal coordinates we can express Rte and Rgg in terms of 
uniformly bounded quantities, as long as R stays bounded away from zero. 



Therefore these functions are uniformly bounded and equation (22) then 
implies that Rtt is uniformly bounded as well. 

In the vacuum case one can take the derivative of the evolution equa- 
tions and repeat the argument in step 2 to obtain bounds on second order 
derivatives of U and A. Here we need another argument. First we write the 
evolution equations for U and A in the forms 

Uu-Uee = { -%^{Ue + U t )-^±^{U t -U e ) 

+ 2& {At " Ae){At + Ae) + 2 6 (56) 

and 

Au-Aee = { -^^(A e + A t) + ^±^l(A t -A 6 ) 

-2{A t - A e ){U + U t ) - 2(A e + A t )(U t - U ) 
+2 Re 2 (n-w)s 23 , (57) 

where k denotes p — P\ + P 2 — P3 ■ Taking the ^-derivative of these equations 
gives 

d x d 6 U e = L +^ U o + ^d x U e 

+^(AxdiA e + A^xAe) + -e 2 ^K 9 , (58) 

and 

dxd^Ae = L + ^A 6 -^dxA e + 2U^dxA 6 + 2A x d^U e 

+2U x d^A e + 2A(.d x U e + 2Re 2 ^~ 2U ) (S 23 ) e , (59) 

Here, L contains only k and S23, first order derivatives of U, A and ??, and 
first and second order derivatives of R, which all are known to be bounded. 
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These equations can of course also be written in a form where the left hand 
sides read d^d\Ug and d^d\Ag, respectively. By integrating these equations 
along null paths to the past of the initial surface, we get from a Gronwall 
argument a bound on 

SUp(|^^| + \d X Ug\ + \d(M + \d X Ag\), 

ees 1 

as long as R is bounded away from zero, under the hypothesis that the 
integral of the differentiated matter terms Kg and (823)9 can be controlled. In 
order to bound these integrals we make use of a device introduced by Glassey 
and Strauss [GS] for treating the Vlasov-Maxwell equation. It is sufficient 
to show how one of the differentiated matter terms can be boundeded since 
the arguments are similar in all cases. Let us consider the integral appearing 
by integrating ( |5g| ) along the null path defined by d\ which involves pg, 



to , 

^- U) v dgf](s,e - (s - t),v)dvds, (60) 



where t G (t, to]. Next, define 



w = V2d x = d t -dg, s = d t + ^dg. 

Hence, dg and dt can be expressed in terms of W and S by 

9 * = ^^(S-W), (61) 

9t = ^r^ {s+ i w) - (62) 

Now, 

[Wf](s, 9-(s- t),v) = d s [f(s, 9-(s-t),v)}, 
and from the Vlasov equation we get 

[Sf](s,e-(s-t),v) = l-K-V v f](s,e-(s-t),v), 

where it is clear which terms have been denoted by K = (K\, K2, K3). By 
using (|6l|) we can now evaluate the integral above by integrating by parts 
(in s for the P^-term and in v for the S-term), so that the remaining terms 
only involve bounded quantities. Note in particular the ^-integrals are eas- 
ily controlled in view of the uniform bound on Q(t). Thus, the integrals 
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of the differentiated matter terms can be controlled and the Gronwall ar- 
gument referred to above goes through. So we obtain uniform bounds on 
\d^Ug\, \d\U$\, \d^Ag\, and |<9,\^0|, and therefore also on \Ugo\, \Ute\, \Agg\ and 
\Ate\, as long as R is bounded away from zero. The evolution equations 



( pOD and (21) then give uniform bounds on \U tt \ and \A tt \. By differentiat- 
ing equation (|23|), it is now straightforward to obtain bounds on the second 
order derivatives of r), using similar arguments to those already discussed 
here, in particular the integrals involving matter quantities can be treated 
as above. Bounds on the first order derivatives of the distribution function 
/ may now be obtained from the known bounds on the field components 
from the formula 

f(t, e, v) = / (6(o, t, e, v), v(o, t, 0, v)), (63) 



since /o is smooth and since <96 and dV (here d denotes dt,dg or d v ) can 
be controlled by a Gronwall argument in view of the characteristic system. 

Step 5. (Bounds on higher order derivatives and completion of the proof.) 
It is clear that the method described above can be continued for obtaining 
bounds on higher derivatives as well. Hence, we have uniform bounds on 
the functions R, U, A, r/ and / and all their derivatives on the interval (t, to] 
if R > e > 0. This implies that the solution extends to t — > — oo as long as R 
stays bounded away from zero. In view of the discussion after the statement 
of Theorem 1, this completes the proof of Theorem 1 in the contracting 
direction. 

□ 



5 Analysis in the expanding direction 

To begin the analysis in the expanding direction (increasing R) in areal 
coordinates we need to start with data on a R =constant Cauchy surface 
(recall that in areal coordinates R = t). That this can be done follows 
from the geometrical arguments in [BCIM] (cf. the discussion following the 
statement of Theorem 1). There it is shown that if Gowdy symmetric (or 
more generally T 2 symmetric) data is given on T 3 , and if Rq is the past 
limit of R along past inextendible paths in D~ on j and if R\ := inf^3 R, then 
for every d G (Rq, Ri), the R = d level set is a Cauchy surface, and these 
Erf foliate the region D~ on j n/ _ (S/j 1 ). Here I~(S) is the chronological past 
of S (see [HE]). The surfaces lie to the past of the initial surface. Let us 
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pick one of them, say S^. The spacetime D~(h,k, fo) induces initial data 
for the areal component fields (U, A, rj, a) and the distribution function / 
on Yit 2= d 2 . By combining the local existence proof in harmonic coordinates 
[CB] , and the arguments in [CI] which show that the spacetime admits areal 
coordinates, we obtain local existence for the initial value problem in these 
coordinates. Now, in order to extend local existence to global existence in 
these coordinates, it is again sufficient to obtain uniform bounds on the field 
components and the distribution function and all their derivatives on a finite 
time interval \p2,tz) on which the local solution exists. 

Step 1. (Bounds on a, U,A and fj.) 
In this step we first show an "energy" monotonicity lemma and then we 
show how this result leads to bounds on fj := n + In a/2 and on U and A. 
Let E(t) be defined by 

E(t)= / [a-mf + ^U^ + ^-,(a-^A 2 t +^A 2 e ) + ^e 2 ^- u ^p]de. 
Js 1 & 



Lemma 1 E(t) is a monotonically decreasing function in t, and satisfies 
j t E(t) = -~ t Jja-^U 2 + ^Aj + ^e^- u \ P + P 3 )]de<0. (64) 



Proof. This is a straightforward computation, using the evolution equation 
for U and A, the constraint equations and the Vlasov equation. First, using 
the evolution equations for U and A we obtain after a short computation 



I / J°-*u? + yftut + ^r(«- 1/2 A 2 + ^Aj)]de 

J S 

+ / ^Ute^-V (P-P1 + P2- P3) + —V^A t S 23 d9. (65) 
Js 1 t 

Next we have 

+ / 2^( Vt - U t )e 2 ^pd9 + / VZe^-WptM. (67) 
Js 1 Js 1 
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Now we use the Vlasov equation and the definition of p to write the last 
integral as 



+ 



s 1 Js 1 

V^e 2 ^-u) \( m -U e + ^)V^(v ) 2 + (Vt - U t )Jv 
L la 



+\faUe[{v 



3\2 l„,2\2 



y/ae 

2U 



] \[olAqv 



2„,3 



t 



df_ 

dv 1 



dvdO 



+ 



+ 



ov z 



nc 



(- - U t )v 3 v° - ^faU B v 1 v 3 ) 



J2U 



+—{v°A t + y/^Agv 



df_ 

dv 3 



dvdO. 



(68) 



Integrating by parts (using the periodicity in 6 and the compact support in 
v) and rearranging the terms we get 



s 1 Js 1 



2{ V -U) 



(v 1 ) 2 . 



,0 



,0 (^) 2 (v 2 ) 2 (« 3 ) 2 x rr e^.uV 
u u v" v u t v u 

3^2 



fdvde 



(69) 



By adding (pq ) and (|67|), using (6£) for the last term in (p7|), we see that 
the terms involving Ut and At vanish and we obtain 



d 
dt 



Eit) 



si , o , /2 V 2 t + «C/, 2 + ^(A 2 + aA 2 ) + aJb-Wp] d6 



Js 1 

f 2U? e 4l/ y/Ei 



e -v, ^( p + p 3 )^. (70) 



By substituting a t and rjt using the constraint equations 



Uf + aU z e +^{Ai + aAj) + e 



2(7,-17) 



ap, 



(71) 
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and 

a t = 2ta 2 e 2( v- u \p 1 - p), (72) 
the first two terms cancel, and ( |64T ) follows. 

□ 

Remark. It is clear from ( |64| ) that a Gronwall argument leads to a bound 
on E(t) also on (0, £2]- F° r T 2 symmetry and vacuum, which is considered 
in [BCIM], this bound is not available. A natural question is then why the 
areal coordinates in our case has to be discarded in the analysis for the past 
direction. However, the analysis of the characteristic system associated to 
the Vlasov equation in lemma 2 depends on the time direction. 



Let us now define the quantity fj by 



fj = r/ + — In a. (73) 
From the constraint equation (Q) we get 

fj e = 2tU t Ue + —A t A e -ty/a~e 2 ( r >- u >J. (74) 

Now, from the elementary inequality \ab\ < ^a 2 + 2cb 2 , for any a,b,c £ 
M, c > 0, and from the fact that \J\ < p, it follows from lemma 1 that for 
any t G [t 2 ,t 3 ) 

/ \fjo\dd <tE(t) <tE{t 2 ). (75) 
Js 1 

Hence, for any 9±,9 2 G S 1 and for any t E [t 2 ,ts) w e have 

\fj(t,9 2 )-fj(t,e 1 )\ 



[ 2 fj d9\ < [ \m\d6 <tE(t 2 ). (76) 
J61 Js 1 



Next, using the constraint equations (Q) and (||), we find that the time 
derivative of fj satisfies 

fj t = t[U 2 + aU 2 + ^{A 2 t + aA 2 ) + ae 2 ^ P 1 ] > 0. (77) 
This relation leads to a control of J sl fjd9 from above, namely 

fj(t, 9)d9 - fj(t 2 , 9)d9 = f t J t V(s, 0)d9^j ds 
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Jt2 Js 1 V« 4s V« 

<sup^a(t 2 ,-) sE(s)ds<C 1 sE{t 2 )ds = CiE(t 2 )(t 2 - t 2 2 )/2. 
s 1 Jt 2 Jt 2 

(78) 

In the first inequality above we used that P\ < p and in the second that a is 
a monotonically decreasing function in t (see (|5|)) and C\ := sup 5 i ^/ a{t 2 , •)• 
We are now in a position to obtain a upper bound on f\ itself. By letting 
C 2 := Jgi fj(t 2 ,6)d6 we get from ( |7"8| ) the inequality 

lc\E(t 2 )(t 2 -t 2 2 ) + c 2 > [ fj(t,e)de (79) 

2 ~ Js 1 



2-7rmaxf7+ / (n-maxn)J. 
Si 7 s i S 1 



(80) 



By applying (^) to the last term we find 

^E(t 2 )(t 2 - t\) + C 2 > 2vrmaxr? - 2irtE{t 2 ). (81) 

Therefore, for some bounded function C(t), we have the upper bound 

maxr7<Cm, (82) 
S 1 

and since fj t > we conclude that r/ is uniformly bounded on S 1 x [£2,^3)- 

Remark. In the analysis below C(t) will always denote a uniformly 
bounded function on ^2^3)- Sometimes we introduce other functions with 
the same property only for the purpose of trying to make some estimates 
become more transparent. 

Next we show that the boundedness of E(t), together with the constraint 
equation (jH|), lead to a bound on \U\. For any 9i, 2 S S l , and t E ^2^3) we 
get by Holder's inequality 

\u(t,e 2 )-u(t,e 1 )\= / u e (t,6)d9 

f f e 2 , \ x /2 , f .e 2 \ 1/2 

< (j a- 1/2 d8j (J ^UjjdOj . (83) 



25 



The second factor on the right hand side is clearly bounded by (Efo)) 1 / 2 ■ 
For the first factor we use the constraint equation (||). This equation can be 
written as 

tya-V 2 ) = t^e 2 ^- u \p - Pi), (84) 

so that for t 6 [£2^3) 



a- 1/2 (t,6) 



sV^e 2(n - u) {p - P 1 )ds + a- 1/2 (t 2 ,e). 



1/2/ 



(85) 



t-2 



Since p> P%, the integrand is positive and bounded by the last term in the 
integrand of E(t). Letting C denote the supremum of a^ 1 / 2 ^,-) over S 1 
we get 



0-2 



a 



V 2 dO < 



t-2 



ae 

s 1 

2 .2 



2( v -u) 



pdOds + 2vrC 



< E(t 2 )(t z -t z 2 )/2 + 2irC. 
Hence, for any 9\,9 2 £ S 1 we have 

\u(t,9 2 )-u(t,e 1 )\<c(t). 



(86) 



(87) 



Next we estimate J sl U(t,9)d9. Let C := j sl U(t2,0)d0, we get by Holder's 
inequality 



U{t,6)d6 
ft 

< 



U t (s,6)d9ds + C 



t 2 Js 



t 2 Js 1 
\U t (s,9)\d9ds + \C\ 



< 



£ Q v 7 ^) 1 (J or x l 2 \J 2 d9^ ' ds + \C\. 



The right hand side is easily seen to be bounded since (|5|) shows that y/a 
is monotonically decreasing and (|64]) gives a bound for the second factor. 
Therefore 



L 



U{t,9)d9 



< C(t), 



for some uniformly bounded function C(t). To obtain a uniform bound 
on U we combine these results. Let U + (t) := maxgi U(t, •), and U-(t) := 
min^i U(t, •). We have 



2nrU±(t) 



U{t,9)d9+ / (U±(t) -U(t,9))d9, 



(89) 
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and the right hand side is bounded from below and above so U is uniformly 
bounded on ^2^3) x S . These arguments also apply to A as well, since the 
factor e 4C/ is controlled by the uniform bound on U. 

Remark. In the case studied in [BCIM], i.e. vacuum and T 2 symmetry, 
a bound on In a, and thus on 77, is directly available. On the other hand, 
the method used here to bound U and A does not apply. This may lead to 
a difficulty in generalizing the result in [BCIM] to matter spacetimes, since 
the bounds on U and A are crucial in order to treat the matter terms when 
the derivatives of U and A are to be bounded. 

Step 2. (Bounds on Ut,Ug, A t , Ag,rj t ,at and Q(t).) 
To bound the derivatives of U we use light-cone estimates in a similar way 
as for the contracting direction. However, the matter terms must be treated 
differently and we need to carry out a careful analysis of the characteristic 
system associated with the Vlasov equation. Let us define 

1 e 4U 
G = -(U? + aUlj) + w (A 2 t +aA 2 e ), (90) 

H = yfaUtUe + jpA t A e , (91) 

and 

d x = -±=(d t + Va~d e ) (92) 

d c = ^Jfr-yfcde) (93) 

A motivation for the introduction of these quantities is based on similar 
arguments as those given in step 2, section 4. For details we refer to [BCIM]. 

Remark. We use the same notations, G and H, as in the contracting 
direction, and below we continue to carry over the notations. The analysis 
in the respective direction is independent so there should be no risk of con- 
fusion. 

By using the evolution equation (0), a short computation shows that 

at 
2^2a 



8 C (G + H) = -^(G + H) 



2rj 

+(U t + V^U e )^=e 2 ^- u \ + ^(A t + Va~A e )S 23 , 
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(94) 



d x (G-H) = *( G -H) 



2y/2a 

AU 



I U? - Va~U t U e + -^(aAj - Va~A t A g ) J 



+([/* - ^)^e 2 ^)« + ^(A t - ^A e )S 23 . 

(95) 

Here k = p—Pi+P2 — Ps. Now we wish to integrate these equations along the 
integral curves of the vector fields d x and respectively (let us henceforth 
call these integral curves null curves, since they are null with respect to the 
two-dimensional "base spacetime" ) . Below we show that the quantity 

r(i) := sup G(t,-) +Q 2 (t), (96) 

is uniformly bounded on ^2^3) by deriving the inequality 

T{t)<C+[ T{s)lnT(s)ds. (97) 
Jt 2 

We begin with two observations. Let 7 and X be a geodesic and a Killing 
vector field respectively in any spacetime. Then 9(7', X) is conserved along 
the geodesic. Here 7' is the tangent vector to 7. In our case we have 
the two Killing vector fields d x and d y . The particles follow the geodesies of 
spacetime with tangent so g^vV^i^xY and g^vT^ipyf are thus conserved. 
Expressing j/ 4 in terms of (see (|i~4|)) we find that 

V 2 {t)e u(t,e(t)) 

and 

V 2 (t)Ae u + V 3 (t)te- U & e ®\ 

are conserved. Here V 2 (t), V 3 (t) and 0(t) are solutions to the characteristic 
system associated to the Vlasov equation. From step 2 we have that U and 
A are uniformly bounded on [£2^3)- Hence |^ 2 (i)| and |^ 3 (t)| are both 
uniformly bounded on [£2^3)1 and since the initial distribution function /o 
has compact support we conclude that 

sup{|u 2 | + \v 3 \ : 3{s,6) e [t 2 ,t] x S 1 with f(s,6,v) ^ 0}, (98) 
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is uniformly bounded on ^2^3)- Therefore, in order to control Q(t) it is 
sufficient to control 

Q 1 ^) := sup-H?; 1 ! : 3(s, 9) £ [t 2 , t] x Shuch that/(s, 9, v) / 0}. (99) 

Below we introduce the uniformly bounded function j(t) to denote estimates 
regarding the variables v 2 and i> 3 . Next we observe that there is some 
cancellation to take advantage of in the matter term (p — Pi) which appears 
in the equations for G + H and G — H above. This term can be estimated 
as follows 

0<(p-Pi)M) = / (v° -^f)f(t,9,v)dv 

1 + (v 2 ) 2 + (v 3 ) 2 - 
/(*> v )dv 

2\2 , ,3\2i' dv 



< [1 + {v ^ + {v -y ] \ f \ 



< II/oIIootW 



dv 1 



'W<QHt) V 1 + ( v1 ) 2 
< C7 7 (t)lnQ 1 (t). (100) 

In a similar fashion we can estimate P2, P3 and 523. Indeed, for k = 1, 2, we 
have 



< Il/ol|oo7(*) 



< C 7 (£) In Q x (t). (101) 

The argument is almost identical for £23. 

Remark. Since the matter of interest is large momenta we have here 
assumed that Q 1 ^) > 2 to avoid the introduction of some immaterial con- 
stants in the estimates. 



Let us now derive (|97|). As in step 2 in section 4 we integrate the equa- 
tions above for G + H and G — H along null paths. For t > £2, let 



A(t,9) = f y/a(s,0)ds, 
hi 
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and integrate along the two null paths defined by d x and d^, starting at 
(t2,0) and add the results. We get for t G [^2^3), 

G(t, 8) = ~[G + H](t 2 , 9 - (A(t) - t 2 )) + \[G + H](t 2 , + (A(t) - t 2 )) 



1 r 


9-(A(s) 


1 /"* 


?-(A(s) 




? - (A(s) 




- (Ms) 



(102) 



where 



Kl = -^-(G + H), K 2 = -^-(G-H), (103) 
1 / _ e iU _ \ 

Ll = "7!t V* + + + v^M*)J » ( 104 ) 

1 / - e 4C/ _ \ 

L2 = "Tit v ' ~ y ^ UtU ° + ~ ^ AtAe n ' ( 105 ) 

M = ^e 2 ^~ u) K, M = e 2i ^S 23 . (106) 

Note that in the expression for M and M we used ae 2r] = e 2f} . It is easy 
to see that both G + H and G — H can be written as sums of two squares. 
From the constraint equation equation (JsJ) we find that at /a < so that K\ 
and K 2 are nonpositive. Using the elementary inequality 2ab < a 2 + b 2 and 
the fact that \fj\ and \U\ are uniformly bounded we obtain from ( 102 ) the 
inequality 

f* 1 

supG(i, •) < supG(t2, •) + supil(i2, ■) + C / - sup G(s,-)ds 
e e e Jt 2 s 

+ f C(s) sup[^/G(s~)((p -P x +P 2 - P 3 ) + S 23 )]ds 
Jt-2 e 

< C + C{t) [ [sup G(s, •) + sup y/G(s, •) lnQ l (s)]ds, 
Jt 2 e e 

(107) 
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where ( [100p and ( |101| ) were used in the last inequality. 

Remark. The sign of K\ and K 2 simplified the estimate above. This 
is not crucial since (c^l/ce is bounded by lnQ 1 (t) which is sufficient for ob- 
taining a bound on Tit). 

Let us now derive an estimate for Q 1 in terms of sup 9 G. 
Lemma 2 Let Q l {t) and G(t,9) be as above. Then 

\QHt)\ 2 <C + D(t) /W(s)) 2 + supG(s, -)}ds, (108) 
Jt 2 

where C is a constant and D(t) is a uniformly bounded function on [t2,ts). 

Proof. The characteristic equation for V 1 associated to the Vlasov equation 
reads 

^ = -to-u e + ?)VSv°-{m-u t )v* 

ds 2a 

-4^((^) 2 " (O 2 ) + ^e 2 W. (109) 
V u sv u 

We will now split the right hand side into three terms to be analyzed sepa- 
rately. Expressing r)g and Tft by using the constraint equations (|3|) and (|j) 
we obtain 

^{V\s)f = 2V\s)^V\s) = T 1+ T 2 + T 3 , (110) 
ds as 

where 

Ti = -2V 1 {s)[sa^- u \jV° + pV 1 )}, 
T 2 = -2V 1 (s) 



e 4C/ 

s (U* + aU% + j^{A$ + aA 2 e ))V l 



+2sy/aU e UtV° - y/a~U 9 V° - U t V l + —^A t A e V° 

2s 



T 3 = -2V\s)l^((V3f-(Vr)-^e 2U V>V\ 
Let us first estimate T\. We split it into two terms 

T\ = Tf + T+ = -2sV 1 {s)e 2 ^- u) {r + J+), (111) 
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where 

I 



= f [ (vW + v°V 1 )f(s,9,v)dv 1 dv 2 dv 3 , 

p poo 

1+ = (vW + v°V 1 )f{s,e,v)dv 1 dv 2 dv 3 . 

Jm. 2 Jo 

Let us now consider the two cases V 1 ^) > and V 1 (s) < 0. On a time inter- 
val where V 1 (s) > 0, I + is nonnegative and T x + can therefore be discarded 
since it is nonpositive. The kernel in Il~ can be estimated as follows 



^ v^-ifiV 1 

(^) 2 (1 + (V 2 ) 2 + (V 3 ) 2 ) {V l ) 2 (l + (v 2 ) 2 + {v 3 ) 2 ) 

Of course, the cancellation of the terms (v 1 ) 2 ^ 1 ) 2 is essential in this com- 
puation. The second term is positive since V 1 (s) > and v 1 < 0, and 
contributes negatively to T-f and can be discarded. The first term is nega- 
tive and the modulus can be estimated by 

(„1)2 (1 + {V 2)2 ± (y3)2) lv l l{1 ± {V 2 ? ± (T/3)2) 

In the expression for T\ we first note that 2sae 2 ( r, ~ u "> = 2se 2 ^~ u ^ < C(s). 
Hence, on the time interval where V 1 (s) > we can estimate T\ by 

< ll/o||C( S ) 7 (s) / t>W<C(s)(QW- (113) 

J 

On a time interval where < we see that is nonpositive and can be 
discarded. We can then estimate by using almost identical arguments 
as for T± and we get also on such a time interval, 

Ti < T+ < C(s)(Q\s)) 2 . (114) 

Let us now consider T2. We again study the cases V 1 (s) > and V 1 (s) < 0. 
Assume first that V 1 (s) > on some time interval. The expression for T2 
can be written T2 = Tf + (p=principal, r=rest) where 

7f = -2(F 1 ( S )) 2 (j s (l/ t + ^) 2 - (C/t + v^)] + [^(^ + V^A e ) 2 }) 
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and 

_ - e 4U _ 

T r 2 = 2(V°(s) - V 1 {s))V 1 (s)[^U e - 2s^U t U e - —^A t A e }. 

For we have 

ini _ ^w tt wr(.) flPri<Wr ^ 

< (s + l) 7 (a) sup G(s,-)- (H5) 



Since the matter of interest is large G we have here assumed that VG< G. 
This assumption will be used below without comment. To estimate T% we 
observe that for s > t%, 

sa 2 — a > — > — , for any a G R. 

The term involving A contributes negatively and can be discarded, thus 

2f < < C{Q\s)f. (116) 

On a time interval where V rl (s) < 0, the same estimates hold. Indeed, we 
only have to write T 2 = + in the form 

T* = -2(V\s)f {[s(U t - V^Ue) 2 - (U t - ^cU e )\ + —(A t - yfcAg) 2 ) 
and 

_ - e 4U - 

T r 2 = 2(V°(s) + V 1 (s))V 1 (s)[\J~aUe - 2s^U t U e - —^A t A e ], 

and the same arguments apply. Therefore we have obtained 

T 2 <T2° + |T 2 r | < C (Q 1 (s)) 2 + C '(s) sup G(s,-). (117) 

e 

Finally we estimate T3. It follows immediately that 



\V\s) 



\T 3 \ < 1 ( s yJ-£>l y fa\Uo + — A e \ <C(s)supG(s,-). (118) 

V S Q 

The lemma now follows by adding the estimates for Tjt, k = 1, 2, 3. 



2U 
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□ 



Combining the estimate for (Q 1 (t)) 2 in the lemma and the estimate (107) for 



sup e G(t, •), we find that T(t) satisfies the estimate (97) and is thus uniformly 
bounded. The constraint equation (|3|) now immediately shows that \rjt\ is 
bounded by 

2tG + te 2 ^- u ^p < C(i)[sup G(t, •) + (Q(t))% 



since 



P= / fdv< H/olloo / dv<C(Q(t)) 6 . 

JR 3 J\v\<Q(t) 

Analogous arguments show that \at\ is uniformly bounded. The uniform 
bound on G provides bounds on \U%\ and \A t \, but to conclude that \Ug\ and 
\Ag\ are bounded we have to show that a stays uniformly bounded away 
from zero. Equation (|||) is easily solved, 

a (t,0) =a(t 2 ,0)e Jtt 2 F{s ' e)ds , (119) 

where 

F(t,9) :=-2tB 2 ^- u \ P -P x ), 

which is uniformly bounded from below. Hence \Ug\ and \Ag\ are bounded 
and step 2 is complete. 

Step 3. (Bounds on df, oeg and rjg.) 
The main goal in this step is to show that the first derivatives of the dis- 
tribution function are bounded. In view of the bound on Q(t) we then also 
obtain bounds on the first derivatives of the matter terms p, J, S23 and Pk, 
k = 1, 2, 3. Such bounds almost immediately lead to bounds on ctg and r]g. 

Recall that the solution / can be written in the form 

f(t, 9, v) = /o(6(0, t, 9, v), V(0, t, 9, «)), (120) 
where 6(s, t, 9,v), V(s, t, 9, v) is the solution to the characteristic system 
d& ^V 1 , . 

^ = -(W - U„ + ^-)V^V° - (r,, - UJV 1 
ds la 

(V 3 ) 2 - (V 2 ) 2 e 2U V 2 V 3 
-V^Ug [V > y0 {V > + G —^Ag^, (122) 
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dV 2 ~ V X V 2 

__ = _ UtV 2 - ^Ug^-, (123) 

dv 3 ,1 „ w , 3 ^ r yV 3 



(- - [Z^ 3 + 



c 



2U 



V 1 . 



^ (A t + V^A e ^)V 2 , (124) 

with the property @(t, t, 9, v) = 9, V(t, t, 9, v) = v. Hence, in order to estab- 
lish bounds on the first derivatives of / it is sufficient to bound dQ and dV 
since /o is smooth. Here d denotes the first order derivative with respect to 
t, 9 or v. Evolution equations for <9G and dV are provided by the character- 
istic system above. However, the right hand sides will contain second order 
derivatives of the field components, but so far we have only obtained bounds 
on the first order derivatives (except for rjg,a$). Yet, certain combinations 
of second order derivatives can be controlled. Behind this observation lies 
a geometrical idea which plays a fundamental role in general relativity. An 
important property of curvature is its control over the relative behaviour of 
nearby geodesies. Let j(u, A) be a two-parameter family of geodesies, i.e. 
for each fixed A, the curve u t— > j(u, A) is a geodesic. Define the variation 
vector field Y := jx(u, 0). This vector field satisfies the geodesic deviation 
equation (or Jacobi equation) (see eg. [HE]) 

D 2 Y 

77T = (125) 

where D/Du is the covariant derivative, R the Riemann curvature tensor, 
and 7' := 7«(n, 0). Now, the Einstein tensor is closely related to the cur- 
vature tensor and since the Einstein tensor is proportional to the energy 
momentum tensor which we can control from step 2, it is meaningful, in 



view of (125) (with Y = dQ), to look for linear combinations of dQ and 
dV which satisfy an equation with bounded coefficients. More precisely, we 
want to substitute the twice differentiated field components which appear 
by taking the derivative of the characteristic system by using the Einstein 
equations. The geodesic deviation equation has previously played an im- 
portant role in studies of the Einstein- Vlasov system ([RR], [Rn] and [R13]). 



Lemma 3 Let Q(s) = Q(s,t,9,v) and V k {s) = V k {s,t,9,v), k = 1,2,3 be 
a solution to the characteristic system ^12\ )- ^12^) . Let d denote dt,dg or 
d v , and define 

$ = a ~ 1/2 dQ, (126) 
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z 1 



dv 1 + 



n t V U t V° (V ) 2 - (V 1 ) 2 + (V 2 ) 2 - (V 3 ) 



a ^ (V ) 2 - (V 1 ) 2 

yl((y2)2_(y3)2 ) At£ 2U V V 2y3 

(V ) 2 - (V 1 ) 2 

v 1 v 2 v 3 



y/at (V ) 2 - (V 1 ) 2 



(y°) 2 - (v 1 ) 2 

dV 2 + V 2 U e d&, 



56, 



p 2U 

(V 3 U e V 2 A 6 )dO. 

s 



(127) 
(128) 
(129) 



Then there is a matrix A = {ai m }, I, m = 0, 1, 2, 3, such that 

n := (tf , Z 1 , Z 2 , Z 3 ) T 

satisfies 

An, 



ds 



(130) 



and the matrix elements ai m = a/ m (s, 6(s), V k {s)) are all uniformly bounded 
on [t 2 ,t 3 ). 

Sketch of proof. Once the ansatz ( |126| )- (|129| ) has been found this is only a 
lengthy calculation. To illustrate the type of calculations involved we show 
the easiest case, i.e. the Z 2 term. 

dZ 2 



ds 



d_ 
ds 



(dV 2 + V 2 U dO) 



= d{^V 2 ) + d ^U e dQ 
ds ds 

+V 2 {U te + Uoe^)d@ + V 2 U e d&. 

ds ds 

(131) 

Now we use ( [1211) and (|123D to substitute for dQ/ds and dV 2 /ds. We find 
that the right hand side equals 



_ v l v 2 ~ v l v 2 

d(-U t V 2 - ^hlUe-^) + (-U t V 2 - y/EUt 



)U e d@ 



+V 2 (U W + Uogyfa^de + V 2 U e (^^o^9 + v^(-j^)) . 

Taking the d derivative of the first term we find that all terms of second 
order derivatives and terms containing ag cancel. Next, since 



_ ■yly2\ /yl 

^ U e d ( -yo~ ) + V^U e V 2 d 



-V^Ue^dV 2 , (132) 



36 



we are left with 



dZ 2 ~ V X V 2 — V 1 

— = -(u t v 2 + yfau e ^-)u e de - (u t + ^u e ^)dv 2 . (133) 

Finally we express this in terms of ^f, Z , Z 2 and Z 3 . Here this is easy and 
we immediately get 

Clearly, the map (8Q,dV k ) \—* (^,Z k ) is invertible so that this step is 
easy also in the other cases. It follows that the matrix elements a2m, rn = 
0,1,2,3, are uniformly bounded on ^2,^3) (only 022 is nonzero here). The 
computations for the other terms are similar. For the Z 1 term we point 
out that the evolution equations (|7j) and (|8|) should be invoked and that the 
matrix element aio contains rj$ and ag/2a, but they combine and form fj. 

□ 

From the lemma it now immediately follows that |f2| is uniformly bounded on 



[t2,ts). Moreover, since the system ( 126 )- (|129| ) is invertible with uniformly 
bounded coefficients we also have uniform bounds on |<90| and k = 

1, 2, 3. In view of the discussion at the beginning of this section we see that 
the distribution function / and the matter quantities p, J, S23 and Pk , are 
all uniformly C 1 bounded. From the constraint equation (||) we now obtain 
a uniform bound on ag by a simple Gronwall argument using as usual the 
identity ae 2<yT, ~ u ^ = e 2 ^ - ^. Finally this yields a uniform bound on 770 since 

and a stays uniformly bounded away from zero. 

Step 4- (Bounds on second and higher order derivatives.) 
It is now easy to obtain bounds on second order derivatives on U and A by 
using light cone arguments. We define G and H by 



1 e 4tv 
G = -(U 2 +aU 2 d ) + w (A 2 t + aA 2 9 ), (134) 

H = V^U tt U t e + -jpAuAto, (135) 

and use the differentiated (with respect to t) evolution equations for U and 
A to obtain equations similar to (|94]) and (|95"|). In this case a straightforward 
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light cone argument applies since we have control of the differentiated matter 
terms. Ugg and Agg are then uniformly bounded in view of the evolution 
equations ([?]) and (|8|). Bounds on second order derivatives on / then follows 
from ( |130| ) by studying the equation for dQ. The only thing to notice is 
that fjgg is controlled by (H). It is clear that this reasoning can be continued 
to give uniform bounds on [tg, £3) for higher order derivatives as well. In 
view of the discussion after the statement of Theorem 1 in section 3, this 
completes the proof of Theorem 1 in the expanding direction. 

□ 
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